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Abstract 

We present a study of W^+multijet events that compares the kinematics of the 
observed events with expectations from direct QCD W+jet production and from pro- 
duction and decay of top quark pairs. The data were collected in the 1992-93 run with 
the Collider Detector at Fermilab (CDF) from 19.3 pb _1 of proton-antiproton collisions 
at y/s = 1.8 TeV. A W+ > 2 jet sample and a W+ > 3 jet sample are selected with 
the requirement that at least the two or three jets have energy transverse with respect 
to the beam axis in excess of 20 GeV. The jet energy distributions for the W+ > 2 jet 
sample agree well with the predictions of direct QCD W production. From the W+ > 3 
jet events, a "signal sample" with an improved ratio of tt to QCD produced W events is 
selected by requiring each jet to be emitted centrally in the event center of mass frame. 
This sample contains 14 events with unusually hard jet E T distributions not well de- 
scribed by expectations for jets from direct QCD W production and other background 
processes. Using expected jet distributions, a relative likelihood is defined and used 
to determine if an event is more consistent with the decay of tt pairs, with M top = 170 
GeV/c 2 , than with direct QCD W production. Eight of the 14 signal sample events 
are found to be more consistent with top than direct QCD W production, while only 
1.7 such top-like events are expected in the absence of tt. The probability that the ob- 
servation is due to an upward fluctuation of the number of background events is found 
to be 0.8%. The robustness of the result was tested by varying the cuts defining the 
signal sample, and the largest probability for such a fluctuation found was 1.9% . Good 
agreement in the jet spectra is obtained if jet production from tt pair decays is included. 
For those events kinematically more consistent with tt we find evidence for a 6-quark 
content in their jets to the extent expected from top decay, and larger than expected 
for background processes. For events with four or more jets, the discrepancy with the 
predicted jet energy distributions from direct QCD W production, and the associated 
excess of 6-quark content is more pronounced. 

PACS numbers: 14.80.Dq, 13.85.Qk, 13.85.Ni 
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1 Introduction 



Recently CDF presented evidence for top quark pair production, both via the observation 
of events with two high Pt leptons and via the observation of events with a W, 3 or more jets, 
and a jet tagged as a b quark jy. In that analysis the distributions of specific kinematic 
parameters of the events, such as jet energies and angles, were not used to discriminate 
between signal and background. It is of interest to search for evidence of top quark pair 
production based on this event structure and to determine whether one can select a top 
quark enriched sample of events with suitable cuts on kinematic variables. 

The main background to ti production comes from higher order QCD production 
of quarks and gluons in association with direct W production. Recent experiments have 
indicated that the top quark mass is larger than of order 130 GeV/c 2 Q g ||. For such a 
heavy top quark it is difficult to distinguish the signal from the background based on the 
properties of the W. However, the jets in ti events have higher energies on average than 
those accompanying the W in direct W production, and are expected to be produced at 
larger angles relative to the beam direction. In this paper we separate ti events where one 
of the VF's decays leptonically and the other hadronically (ti — ► W + b + W~b — > Ivb + qqb) 
by exploiting these properties. 

The paper is organized as follows. Section 2 describes the analysis cuts used to define 
the W+ jet sample. Section 3 gives a brief comparison of the kinematics of directly produced 
W+ jet events and of W+ jet events from top quark decay, and explains the analysis 
strategy. Section 4 summarizes various comparisons of QCD Monte Carlo predictions which 
successfully fit experimental measurements in processes where a top quark contribution 
can be neglected. Kinematic features of our W+ > 3 jet data sample are compared to 
background and to top quark prediction in Section 5. This comparison shows evidence for a 
top quark-like component in the data. Section 6 combines this result with the independent 
information obtained from the algorithms which provide identification of b quarks in the 
events. In Section 7, as an additional test, we look for an excess of W + 4 jet events in the 
top quark candidate sample. The conclusions are presented in Section 8. 
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2 Data Selection 

This analysis is based on 19.3 pb -1 of data from 1.8 TeV pp collisions taken with the 
CDF detector during the 1992-1993 Tevatron run. The CDF detector is described in detail 
elsewhere Q , 0j , |§ • For this run, the tracking system was upgraded with a high precision 
silicon vertex detector (SVX) ||), and the muon detector was improved at pseudorapid- 
ity @ M<0.6 by adding an absorber of 0.6 m of steel followed by drift chambers. In 
addition, the coverage of the central muon detector was extended to the region of pseudora- 
pidity 0.6 < |?7| < 1.0 (over about 2/3 of the azimuth) with drift chambers and scintillation 
counters. The transverse momentum P^, defined as = Psin#, is the projection of the ob- 
served momentum (P) onto the plane transverse to the beam axis. Similarly, the transverse 
energy is defined as: E^= Esin#, where E is the energy measured in the calorimeter. The 
identification and measurement of isolated, high-Py electrons and muons, the measurement 
of the missing Ej- ($ T ) indicative of neutrinos in the events, the jet clustering algorithm, 
and the jet energy corrections are discussed in Ref. |l|,0 and (8|. 

A sample of W — > ev{[iv) candidate events was selected with the requirement that 
Ef,>20 GeV (P^>20 GeV/c) and $ T >25 GeV. In addition, the transverse mass, defined 
as M.t = pEj-Fj-rU-cosAc/))] 1 / 2 (where A(j) is the difference in azimuthal angle between the 
missing energy direction and the lepton), was required to be larger than 40 GeV/c 2 . The 
jets are reconstructed with a cone size R = \J A<I> 2 + At/ 2 = 0.4 (where A$ is the cone 
half-width in azimuth and Ar/ is the cone half-width in pseudorapidity) . The jet energies 
are corrected by a rapidity and energy dependent factor which accounts for calorimeter non- 
linearity and reduced response at detector boundaries ||. In addition to these detector 
effects, a correction is also made for energy which is radiated out of the jet reconstruction 
cone. The % T is calculated after correcting the jet energies. In order to allow for a clean 
separation of jets from each other and to facilitate the comparison of energy distributions 
with theoretical expectations, the centroids of the three leading jets are required to be sep- 
arated from each other by AR > 0.7. 

Backgrounds from Z decay, Drell Yan production of dileptons, and possible ti events 
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in which both W J s decay leptonically are removed by rejecting events with an additional 
isolated track with P*r > 15 GeV/c in the central tracking system that is not associated with 
the primary lepton. Tracks are defined as isolated when the total transverse momentum of 
the charged tracks (other than the track in question) in a cone of radius R = 0.4 centered 
on the track is less than 0.1 times the Pt of the track. A study of a QCD multijet sample 
has shown that fewer than 1% of jets with E^(jet) > 20 GeV are rejected by this cut Q. 
An additional Z removal algorithm eliminates events with an oppositely-charged dilepton 
(ee or fi/j,) invariant mass in the range 70 to 110 GeV/c 2 . 

A sample of W+ > 2 jets with the two leading jets having EyQet) > 20 GeV and 
|?](jet)| < 2 is selected, where the tt contribution is expected to be relatively small. This 
sample is studied in order to check whether the energy spectra of the leading jets agree with 
QCD prediction for direct W+ jet production. The search for a tt component is performed 
in the subsample with > 3 such jets. 

The primary differences in event selection between this analysis and the analysis of 
W+ > 3 jets performed in Ref. [Q| are that (1) corrections to jet energy and Ex are made 
prior to event selection, (2) the jets are explicitly required to be separated by Ai? > 0.7, 
(3) a cut is added on the transverse mass, and (4) the rejection of events with an additional 
isolated track is included. These changes are made to simplify the comparison of observed 
jet energies with theoretical predictions for direct W+ jet production and to reduce back- 
ground from non-VK sources. 

The fraction of all ti events that should fall into our W+ > 3 jet sample is deter- 
mined from the HERWIG pp — ► tt Monte Carlo program and the CDF detector simulation. 
Corrections to the acceptance for trigger inefficiencies and differences in lepton identifica- 
tion between data and Monte Carlo simulation are identical to those described in Ref. |l[. 
For this analysis we find that the top quark acceptance ranges from 2.7±0.2% at Mt op =150 
GeV/c 2 to 3.0±0.2% at M top =190 GeV/c 2 . The number of tt events expected in the W+ > 3 
jet sample using the Standard Model top quark production cross section from Ref. [|ll]] is 
about 7 for M top = 170 GeV/c 2 . Using the cross section from Ref. [[jj we expect 16.71 6 'q 
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events. We observe a total of 49 events, 25 of them being in common with the 52 W+ > 3 
jet event sample of Ref. [Q]. For top quark events the two sets of cuts make little difference: 
about 90% of all top quark events contained in the sample of 49 events will also show up 
in the sample of 52. QCD W+ jet events often will be close to the jet cuts. Therefore 
only approximately 67% of QCD events from the sample of 49 will also be found under the 
cuts of the sample of 52. If we assume that all 49 events are from QCD, then we expect an 
overlap of approximately 33 events to be compared to the observed overlap of 25 events. 

3 W and ti Kinematics 

Monte Carlo event samples are used to compare the distribution of several kinematic 
variables for top quark and background events. Samples of top quark events of various 
masses were generated with both ISAJET and HERWIG [|l^], and it was verified that 
both Monte Carlo generators give similar results. W+ jet events were generated according 
to the lowest order matrix elements for the production of a W with n final state partons. The 
complete sets of matrix elements at tree level have been determined for n = 0, 1, 2, 3, 4 and 
are implemented in the program VECBOS [fLS] ]. To avoid infrared divergences which would 
occur at small angles and small Pp, cuts are applied in the event generation that require Pt 
(parton) > 10 GeV/c, \rj (parton)| < 3.5, and AR (parton-parton) > 0.4. Unless otherwise 
noted, Q 2 = M^ has been used for the a s scale and the structure functions; this choice 
yields the hardest jet energy spectrum of a number of Q 2 scales considered. Two different 
techniques are used to transform the partons produced by VECBOS into hadrons and jets, 
which can then be processed through the CDF detector simulation. One method employs 
a Field and Feynman fragmentation function |L4|] ("SETPRT"), tuned to reproduce the 
features of observed inclusive QCD jets. The other ("HERPRT" |l5|]) uses part of the QCD 
shower evolution Monte Carlo HERWIG. In this case the events generated by VECBOS 
are assigned an appropriate flavour and colour configuration, and are processed through 
the HERPRT initial and final state evolution program. Unless otherwise noted, the results 
presented here will use the HERWIG approach. The Monte Carlo events have then been 
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processed through a full simulation of the CDF detector and reconstructed in the same 
manner as the data. 

The choice of suitable kinematic parameters to distinguish top quark events from 
background is presently the subject of considerable work. The DO collaboration has used 
event aplanarity and the scalar sum of the jet transverse energies Q. CDF has studied these 
variables as well as other parameters involving combinations of jet energies and angles fl6| . 
Work is presently in progress to identify which parameters provide optimal information. In 
this study we have focused the analysis on jet transverse energies and polar angles. Our 
studies have indicated that these variables are among the most powerful at separating top 
quark signal from direct W+ jet background. 

Jets are ordered in Et with jeti having the highest energy, Er\. It was found that 
the Et2 or Et3 variables are better discriminant between QCD background and top quark 
events than Eti ■ A qualitative indication of the separation that can be obtained between tt 
and direct W+ jet production on the basis of Et2, Et3 is shown in Figure 1, which presents 
the predicted density ±d 2 a/dE T2 dE T3 of W+ > 3 jet and tt events (M top = 170 GeV/c 2 ). 
The distributions are different for heavy top quark and background events, with tt events 
characterized by higher Et jets. 

Selection of events based on the presence and energy of a fourth jet is also predicted 
to be a good discriminant between tt and direct W production. However, in this analysis 
we do not require a fourth jet. This is done in order to minimize: (a) uncertainties in the 
theoretical calculation of the spectrum in tt events with accompanying gluon radiation, 
(b) the uncertainty in the reconstruction efficiency and in the energy measurement of low 
energy jets. The presence of a fourth jet will later be examined in this paper as an indication 
of whether top quark is present. 

Another variable which can discriminate between W+ jets and ti is |cos#*| ||17|| , 
where 9* is the angle between a jet and the incident proton direction in the center of mass 
of the hard subprocess. The component of the hard subprocess center of mass velocity 
along the beam direction is calculated using the four-momenta of the W and all jets with 
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(o) QCD (b) Top 




Figure 1: £d 2 <7 / dE T2 dE T3 for (a) QCD W + 3 jet and (b) top quark (M top = 170 GeV/c 2 ) 
Monte Carlo events. 
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Et > 15 GeV. Jets are included down to this relatively low energy in order to reconstruct 
the laboratory velocity of the initial state subprocess as well as possible. In calculating 
the W 4-momentum, the longitudinal component of the neutrino cannot be determined 
unambiguously and for simplicity is taken to be zero. The expected distribution of the jets 
as a function of | cos9*\ max , the maximum of |cos0*(jetj)|, i=l,2,3 is shown in Figure ||(a). 
The inclusive jet distribution for direct W events is peaked in the forward direction while 
that for top quark events is more central. As in Ref. jeti, jet2 and jet3 are required to 
have |?y(jet)| < 2. The | cos#*| distribution after this cut is shown in Figure §(b). After the 
|r?(jet)| < 2 cut, our studies indicate that a | cos 6*\ cut still improves the signal/background 
ratio. It also allows one to define a background depleted "signal sample" as those events 
in which each of the three leading jets satisfies |cos#*| < 0.7, and a background enriched 
"control sample" which contains all events in which at least one of the jets has | cos 6*\ > 0.7. 
The Monte Carlo predictions show that the |cos 9*\ cut generates a harder jet distribution 
for top quark production, while for direct W+ jet production it leaves the Ey distributions 
essentially unaffected. Therefore an analysis which attempts to separate top quark from 
background based on the shape of the Ey distributions can be expected to become more 
discriminating after applying the |cos#*| cut. 

Using the Standard Model tt cross section from Ref. ]l(| we expect approximately 6 
events in the signal sample and 7 in the control sample for a top quark mass of 150 GeV/c 2 , 
while for a top quark mass of 190 GeV/c 2 we expect approximately two events each in the 
signal and control samples. Assuming the top quark production cross section from Ref. Q, 
the number of top quark events expected for Mt op = 170 GeV/c 2 is 7.7^2'g i n the signal 
sample and 9.0l 3 ' 2 in the control sample. The signal sample contains 14 events and the 
control sample 35 events. Therefore in the signal sample a signal/background of the order 
of 1 can be expected, while in the control sample this ratio would be nearly three times 
worse. 
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Figure 2: Distributions of the | cos 0*\ max variable predicted by the HERWIG top quark 
(M top = 170 GeV/c 2 ) and VECBOS 1T + 3 jet calculations for: (a) the inclusive distribution 
and (b) after applying the cut on |?7(jets)| < 2. The distributions are normalized to unit 
area. 
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4 Reliability of W+ Multijet Predictions 

In subsequent sections a detailed comparison is made for the observed jet energy dis- 
tributions for events that contain a W and > 3 jets with the predictions of QCD direct 
production of W and jets (as implemented in the VECBOS program). It is therefore im- 
portant to investigate to what extent these predictions are reliable. Previously CDF has 



compared the cross section for W + n jet production (n < 4) with QCD predictions [18] and 
found good agreement. In addition, the jet energy distributions and rapidities for W + 1 jet 
and W + 2 jets show good agreement with the QCD calculations |H| ]. The UA2 collabora- 
tion at CERN has examined the transverse energy distributions for multijet events with up 
to six final state partons and found good agreement [|l9| with expectations from QCD. CDF 
has also found excellent agreement between observation and QCD predictions for inclusive 
distributions in 3 and 4 jet data samples Although the UA2 comparison and the CDF 

multijet comparison involve a different set of matrix elements than for jets associated to W 
production, they demonstrate that in terms of jet detection and reconstruction, excellent 
agreement is obtained between observations and theory for events containing as many as 
four jets. 

A good test is provided by the CDF W+ > 2 jet data sample, which has relatively 
high statistics, is kinematically similar to the W+ > 3 jet sample, and has a relatively small 
fractional contribution from top quark. 

The angular distributions of data and VECBOS events are compared in Figure || in 
terms of the variable | cos0*\ max , which is here the maximum of |cos#*(jetj)|, i=l,2. As for 
the W+ > 3 jet sample, a cut on the jet rapidity was applied at |^(jet)| < 2. The Monte 
Carlo prediction is normalized to the data. The agreement is excellent. The En and E^2 
distributions of these W+ > 2 jet events are compared to VECBOS predictions in Figure |] 
under the requirement that both jets have |cos#*| < 0.7. The Monte Carlo prediction is 
normalized to the data. The confidence level of the likelihood that the data are consistent 
with Monte Carlo predictions is 55% for Fig. 4(a) and 69% for Fig. 4(b). While the 
VECBOS calculation has no phenomenological parameters, the results do depend on the 
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Figure 3: | cos 6*\ max for W+ > 2 jets showing data (points with error bars) and VECBOS 
W + 2 jet events (histogram). 
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Figure 4: distributions of W+ > 2 jet showing data (points with error bars) and 

VECBOS W + 2 jet events (histograms), (a) leading jet, (b) second to leading jet. 



16 



choice of the Q 2 scale and on the minimum separation and Pt of the generated jets. Use 
of Q 2 =< Pt > 2 instead of Q 2 = yields softer spectra. The agreement between data 
and predictions is equally good. With the choice Q 2 =< Pt > 2 the confidence level of the 
likelihood that the data are consistent with Monte Carlo predictions is 62% for both the 
Eti and Et2 distributions. Since top quark events are expected to give harder jet energy 
spectra than direct W+jet production, our default choice, Q 2 = Myy, (see Section 3) is 
conservative. Additional tests on the sensitivity to the Q 2 scale of W+ > 3 jet Monte Carlo 
predictions are shown in the next section. Finally, a test of the predictions of jet production 
associated with vector bosons may be obtained from the Z+ > 3 jets sample, where little 
contamination from Standard Model top quark events is expected. Data and predictions 
are shown in Figure ||. While the statistics are limited, the agreement is good. In the higher 
statistics Z+ > 2 jet sample, one also finds good agreement between data and predictions 
for the Et distributions of the two leading jets. 

5 Kinematic Analysis 

As discussed in Section 3, a signal sample of 14 events is defined by the requirement 
that the three leading jets have |cos#*(jet)| < 0.7. The background enriched control sample, 
where at least one jet has |cos6>*(jet)| > 0.7 contains 35 events. Figure ^ shows the Et 
distributions of the three leading jets in the signal enriched sample. Figure shows the 
same distributions for the process W + 3 jets as predicted by VECBOS and for top quark 
production modeled with the HERWIG Monte Carlo at a top quark mass of 170 GeV/c 2 . 
The distributions are normalized to unit area. The Et distributions of the data are harder 
than those expected from VECBOS. To combine the information from both Et2 and Et3, 
a discriminating function, "absolute likelihood", is defined as follows: 



L a bs — {— 



1 da 



1 da 



(1) 



a dET2 



a dET3 
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Figure 5: distributions of second and third leading jets in Z+ > 3 jet (inclusive) data 
(shown as points with error bars) and VECBOS Z + 3 jet events (histograms). 
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Figure 6: Jet energy distributions for the three leading jets in the 14 events passing the 
signal sample selection cuts. There is one overflow in E^i at Eti = 224 GeV. 
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Figure 7: Jet energy distributions for HERWIG top quark (solid line) and VECBOS 
W + 3 jet events (dashed line) passing the signal sample selection cuts. Each distribution 
is normalized to unit area. 



20 



that is, as the product of the two differential transverse energy distributions each normalized 
to unit area. The are derived from the Monte Carlo simulated distributions fitted by 
analytical functions. A L a b s function can be defined for any process for which a model 
exists, in particular for QCD W+ > 3 jets (L^" ) and top quark (L*^). The L a b s 's can be 
combined to define a "relative likelihood" (L re /) for top quark versus QCD as 

Lrel = 1&/L% D . (2) 



Note that the absolute likelihoods are not probabilities, since E^2 and Et3 are correlated. 
The relative likelihood allows one to compare each individual event to the expectation from 
QCD and from top quark in terms of a single number. This "kinematic tag" provides 
a natural definition of the cut which discriminates events which are more top quark-like 
from events which are more QCD-like. A possible disadvantage of L re ; is its dependence 
on Mtop- In particular, a L re ; which is optimized for a certain top quark mass may have 
reduced sensitivity if the actual top quark mass is significantly different from the assumed 
mass. We choose for our analysis Mt op = 170 GeV/c 2 , based on the results of Ref. Q. We 
discuss the effect of a possible different choice in Section 5.1. 

5.1 Data— Monte Carlo comparison 

Figure 8 shows a comparison of the expected and observed distributions for ln(L^ D ) 
for the W+ > 2 jet and W+ > 3 jet signal and control samples. In the case of W+ > 2 
jets, the L a b s is defined as the product of the E^i and distributions, since a third 
jet is not always present in the event. The W+ > 2 jet sample is expected to have a 
small top quark fraction. The comparison with the VECBOS prediction (Fig. 8(a)) shows 
good agreement. The W+ > 3 jet control sample data, where the QCD background is 
expected to dominate, also agree with the QCD prediction as shown in Figure ||(b). In 
the W+ > 3 jet signal sample (Figure ||(c)), where a ti contribution could be present, 
VECBOS predictions and data are somewhat different. In order to check how significant 
this difference is we performed a likelihood calculation, by assuming a Poisson distribution 
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Figure 8: QCD-predicted absolute likelihood distributions (histograms) compared to data 
(points with error bars), (a): the two leading jets in W+ > 2 jets, a cut |cos#*| < 0.7 was 
applied to jeti and jet2; (b): jet2, jet3 in the control sample; (c) : jet2, jet3 in the signal 
sample. 



with mean equal to the Monte Carlo prediction for each bin. A confidence level of 19.6% is 
found for VECBOS to agree with the data of Figure |(b), and 2.7% for Figure |(c). The 
likelihood for Figure ||(c) is small enough to suggest a component in the signal sample not 
well described by VECBOS. The expectations for HERWIG tt events, when interpreted as 
direct W+ jet events, are shown in Figure || for a number of different top quark masses. 
The signal sample data distribution in \n{L^ s D ) seen in Figure ||(c) is consistent with a 
combination of direct W+ jet events and of tt events in a wide mass range around M( 0p 



= 170 GeV/c 2 . Figure |T|(a) shows how the VECBOS W + 3 jet events and HERWIG 
top quark events are distributed in relative likelihood, ln(L*]J ), when the cuts of the signal 
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Figure 9: Expected distributions of ti Monte Carlo (HERWIG) events as a function of 
ln(L^ D ), when interpreted as direct W+ jet events. The distributions are normalized to 
unit area. 
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sample are applied. The symbol L']J indicates that M top = 170 GeV/c 2 was used to predict 
the expected E^2 and distributions. The two distributions are separated well enough 
to make a top quark signal visible, provided the signal/background is of order 1 as argued 
in Section 3. Figure [lO|(b) shows how the data events of the signal sample are distributed in 
ln(L*gJ ), along with the VECBOS distribution from Figure fK](a), normalized to the data 
at ln(L*gJ°) < 0. The data are not distributed as expected from a pure QCD W+jet sample, 
and look like a superposition of tt and QCD events. 

We find 6 events with ln(L*gJ°) < (more QCD-like) and 8 events with ln(L*gj°) 
> (more top quark-like). If we normalize the VECBOS distribution, 78% of which is 
expected to have ln(L*gJ ) < 0, to the 6 events observed in that region, then we expect 
1. 68^0^2 VECBOS events with ln(L*gJ°) > compared to 8 events observed. This excess 
represents kinematic evidence for the presence of tt production. 

Although M top = 170 GeV/c 2 was assumed, the result is not sensitive to the precise 
value of M top . As an example, Figure |ll|(a) shows the expected distributions in ln(L*gf°) 
for VECBOS W + 3 jets and HERWIG events if M top = 150 GeV/c 2 is assumed. The two 
samples are still well separated. The data are distributed as shown in Figure 0(b), and 
still indicate a superposition of the two processes. The results of a two-component fit to 
the m(L^J°) distribut ion as a function of Mt op and the significance of the observed excess 
at positive ln(L^J°) will be examined in a later section. 

5.2 Evaluation of the statistical significance 

The probability that a background fluctuation can produce the observed excess at 
m(L*g| ) > is calculated from the binomial probability that given the 14 signal sample 
events, they are distributed with 8 or more events in the positive ln(L']J ) side. The cal- 
culation makes use of the fraction of QCD W+ jet events expected at hi(L**J°) < and 
of the statistical error on the fraction. For the primary result (shown in line 1 of Table 1) 
VECBOS W + 3 parton production with a choice of Q 2 = M 2 ^ and the HERPRT fragmen- 
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Figure 10: Ln(L* 1 J°) for QCD VECBOS, top quark HERWIG {M top = 170 GeV/c 2 ) and 
data events for W + 3 or more central jet events (signal enriched sample), (a): W+ > 3 
jet VECBOS (dotted histogram) and top quark HERWIG (solid histogram), normalized 
to unit area, (b): data ( solid histogram) and VECBOS (dotted histogram). VECBOS is 
normalized to data in the region ln(L*gj°) < 0. 
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Figure 11: (a) Expected distribution in ln(L*gf°) of QCD and ti events, for M top = 150 
GeV/c 2 ; (b) the data displayed as a function of the same variable. The dotted histogram 
shows the QCD distribution normalized to the data at ln(L*gf°)<0. 
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tation is used. In this case the probability to observe 8 or more events from a background 
fluctuation is 0.5%. This disagreement between observation and the VECBOS prediction 
is large enough to suggest the possibility that either VECBOS is wrong, or that there is an 
additional high Et process present in the data sample. 

Next we address systematic uncertainties in the Monte Carlo predictions. One sys- 





% rate at ln(L* 1 J u )< 


1. HERPRT (Q 2 =M^) 


78.2±2.4 


2. Et scaled down 


81.2±2.6 


3. Et scaled up 


74.9±2.2 


4. HERPRT (Q 2 =<P T > 2 ) 


80.4±5.5 


5. SETPRT (Q 2 =<P T > 2 ) 


81.1±2.8 


6. HERPRT (Q 2 =M^) 

+ Systematics + non-VF backgrounds 


78±5 



Table 1: Fraction of background events at ln(L*gj°) < for different predictions of the 
ln(L*gj°) shape. The predictions compare to 6 data events observed at ln(L*gJ°) < and 8 
events at ln(L*^°) > 0. 

tematic uncertainty is the possible difference in the energy scale of data and Monte Carlo 
events. The relative uncertainty in the jet energy scale of the calorimeter decreases with 
increasing the jet energy. These effects have been taken into account by varying Monte 
Carlo jet energies by an uncertainty from ±10% at 8 GeV to ±3% at 100 GeV to account 
for detector effects, in quadrature with a ±10% uncertainty due to the assignment of ener- 
gies to partons in the presence of gluon radiation, which is the dominant uncertainty. The 
second and third lines of Table 1 show the results. The uncertainty in VECBOS due to the 
lack of higher order contributions can be addressed by changing the Q 2 scale in a s . For 
comparison with the results shown in the first line of Table 1, the results for Q 2 =<Pt> 2 
are shown in the fourth line of Table 1. The fourth and fifth lines compare the results using 
our default fragmentation algorithm HERPRT with SETPRT (see Section 3) and show very 
little difference. 

Contributions to the event sample from background sources other than the dominant 
direct PF±jet production were studied to determine if they could explain some of the excess 
at ln(L*gJ°) > 0. These additional backgrounds are of two types. First, in the VF±jet sample 
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there is a fraction of non-VK events (e.g.: hadrons misidentified as electrons or muons, or 
real leptons from bb). As in Q, the number of such events is estimated by extrapolating 
the number of events which pass the E T cut but have non-isolated leptons, to the region 
in which lepton isolation is required. When only the isolation cut is released in the signal 
sample, no additional event enters. Following this procedure, the signal sample is estimated 
to contain O.OIq'q events from this source. The ln(L*gJ°) distribution of the non-W^ events is 
shown in Figure |i~2|(a) and is similar to the Monte Carlo predicted distribution of VECBOS 
events. 

A second background is WW, WZ events or single Z events with one non-identified 
decay lepton. The ISAJET Monte Carlo is used to simulate these backgrounds. The WW 
and WZ backgrounds are normalized with the next -to-leading order computations of the 
cross section from Ref. ^0|. For the Z background the normalization is provided by the 
measured CDF Z — > ee cross section [ffijl . The estimated number of such events in the 
signal sample is 0.9±0.3 WW events, 0.13±0.05 WZ and 0.14l^ misidentified Z's. The 
ln(L*g[°) distribution of the dominant WW contribution is shown in Figure |l2|(b). Again, 
most events are at ln(L*JJ°) < 0, as for the QCD single W+jet background. (21±9)% of 
non-W events and (20±4)% of WW events are expected at ln(L*gj°) > 0, compared with 
the (21.8±2.4)% for the W+ jet background (top line in Table 1). 

The probability that the observed excess at positive ln(L*]J ) is consistent with the 
VECBOS prediction including the effects of non- VECBOS backgrounds and the other sys- 
tematic errors discussed above is computed as follows. Non- VECBOS events are chosen 
from a Poisson distribution with the means presented above, and are distributed at positive 
or negative ln(L*gJ°) according to the determined fraction. The remainder of the 14 events 
are taken to have the ln(L*gJ°)<0 fraction predicted by VECBOS, which is taken to be (78 
± 5)%. As can be seen from Table 1 this adequately allows for the variations due to changes 
in the energy scale, Q 2 scale, and the statistical error. The probability is calculated via a 
Monte Carlo program that includes all the uncertainties mentioned above. The probability 
is 0.8% to observe > 8 events with ln(L*gJ°) > in a sample of 14 events originating from 
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Figure 12: (a) Ln(L£J°) distribut ion of non-iy background, as derived from the study 
of a sample of non isolated leptons, with small $ T ; (b) Monte Carlo predicted ln(L*gj°) 
distribution of WW events. The distributions are normalized to unit area. 



direct W+ jets and non-W sources. 

We tested whether the results are stable under reasonable variations in the event 
selection requirements for the signal sample. When we change the requirement on En from 
20 to 50 GeV, change the cut on cos#* from 0.7 to either 0.65 or 0.75 or change the jet-jet 
separation cut from AR = 0.7 to 0.6, we get the probabilities for a statistical fluctuation 
of 0.5%, 1.9%, 0.5% and 0.7% respectively. In the worst case the background fluctuation 
probability is 1.9%. As an additional test, events were selected with the requirement that 
the uncorrected missing transverse energy # T raw > 20 GeV, and no cut on the transverse 
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Figure 13: Distributions of ln(L*gJ°) for the control sample, (a) Distributions of VECBOS 
VK + 3 jet (dotted histogram) and HERWIG top quark events (solid histogram), normalized 
to unit area, (b) 35 data events (solid histogram) versus VECBOS (dotted, with statistical 
errors). VECBOS has been normalized to data in the region ln(L*gJ°) < 0. 

mass of the W . This is the selection used in Ref. This sample has about 3 times larger 
background from fake W events due to misidentified leptons, and the neutrino transverse 
momentum for real W events is not as well determined. However, the acceptance for W 
and top quark events is ~25 % larger. This results in a signal sample of 19 events; 11 events 
are at ln(L£J°) < 0, 8 at ln(L* x J ) > 0. The probability that > 8 events have ln(L^J°) > 
if the events were entirely QCD background is 3.8%. 



Figure 13 shows how the 35 events of the control sample are distributed in ln(L*g] 



together with the predictions from VECBOS and HERWIG. The data are well described 
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by the VECBOS QCD expectation. There is no statistically significant indication for an 
excess at ln(L*]J ) > (whether this is reasonable is addressed in the next section). For 
this comparison, Q 2 =M^ is used as the scale for a s and structure functions. The use of 
Q 2 =<Pt> 2 would give a softer spectrum in ln(L*gJ ) and the number of predicted QCD 
events at ln(L£J°) > would fall by 25%. 

5.3 Cross Section Calculation 

We assume here that the excess of high jet events in the signal sample results from 
ti production and decay. The most probable numbers of ti (Nj op ) and W (Nw) events 
observed in the signal-enriched and background-enriched samples are estimated using an 
unbinned maximum likelihood fit to the observed ln(L*gJ°) distribution. For this fit, we use 
the ln(L*gJ°) shapes predicted from the HERWIG Monte Carlo program for it events (M top 
= 170 GeV/c 2 ), and from the VECBOS Monte Carlo program for QCD W-\- jet events; 
these predicted shapes are shown in Figure |i~(i|(a). The systematic uncertainties on Nt op 
and the ti production cross section, a^, are estimated as a function of Mt op . The following 
effects are considered: (1) uncertainty in the jet energy scale, estimated as described in 
Section 5.2; (2) uncertainty in the Q 2 scale employed by VECBOS to determine the jet E^ 
spectrum in W events^, estimated by choosing Q 2 = < Pj- > 2 rather than Q 2 = M 2 ^; (3) 
Monte Carlo statistics and uncertainty on the lepton detection efficiency; (4) choice of ti 
Monte Carlo generator (HERWIG, PYTHIA, ISAJET); (5) change of the ln(L rei ) shape due 
to variations of assumed top quark mass in the range 150 < Mt op < 190 GeV/c 2 ; (6) the 
inclusion of non-W, WW, WZ, and misidentified Z contributions in the likelihood fit; (7) 
uncertainty in the fitting procedure; and (8) uncertainty in data integrated luminosity (this 
enters only in the calculation of a^) . The results for the signal sample are summarized in 
Table |2| for a^. The systematic uncertainties on Nj op are similar to those on <r^; only the 
totals are listed in Table ||. The number of tt events is found to be N top =6.4l3 2 and 
Nfop = 0.81q '| tog f° r the signal and control samples respectively, where the first error is 

1 We address the effect of the Q 2 scale on the VECBOS ln(L"{°) shape only. Absolute rate predictions 
do not affect the likelihood procedure. 
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150 GeV/c 2 


170 GeV/c 2 


190 GeV/c 2 


(1) Jet Ej scale 


+21% -19% 


+22% -11% 


+23% -5% 


(2) VECBOS Q 2 


+10% -0% 


+10% -0% 


+10% -0% 


(3) MC stat + lepton e 


±10% 


±9% 


±9% 


(4) tt generator 


±10% 


±9% 


±10% 


(5) non-Ty background 


±2% 


±2% 


±2% 


(6) fitting procedure 


±3% 


±3% 


±3% 


(7) Luminosity 


±3.6% 


±3.6% 


±3.6% 


Nfop total 


+27% -21% 


+28% -17% 


+24% -16% 


(jfj total 


+31% -24% 


+33% -23% 


+29% -20% 



Table 2: Individual systematic errors on in the signal sample as a function of M.t op 
are listed in rows (1) to (7). Total systematic uncertainties for both N top and are 
summarized in the final two rows. 

statistical and the second error is systematic. 

The fits indicate more tt candidate events in the signal-enriched sample than in 
the background-enriched sample. The ratio of top quark events in the control sample 
to top quark events in the signal sample predicted by the Monte Carlo calculation is 1.17 
(9.0/7.7). The data fit finds 0.13. However, the statistical significance of the difference from 
expectation, taking into account the errors, is within 1 a. The systematic and statistical 
errors in the determination of Ni 0p are significantly larger in the control sample, due to 
the larger number of QCD P^+jet events. From N top , we calculate the corresponding tt 
cross section. This analysis is performed on the signal sample to minimize the systematic 
effects from the uncertainties in the prediction of the ln(L*]J ) shape for QCD VF+jet 
events. Table || shows the total tt acceptance, including branching ratios, lepton detection 
efficiencies and energy scale uncertainty, as a function of the top quark mass, and the 
results for both Nt op and crrz. Figure [L4| shows the summed ln(L*gJ°) distribution for tt 
and VECBOS corresponding to the Table 3 result (with Mt op = 170 GeV/c 2 ) compared to 
the observed ln(L*gJ°) distribution. The cross section determined from the signal sample 
is consistent with that found in and, given the large errors, is not inconsistent with the 
value of 5.7i ; 6 pb predicted by the theory for a top quark mass of 170 GeV/c 2 [{Tcfl . 
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M top =150 GeV/c 2 M top =170 GeV/c 2 M iop =190 GeV/c 2 


tt acceptance 
N top 


2.7±0.2% 2.9±0.2% 3.0±0.2% 

7 q+4.6 +2.1 fi 4 +3.8 +1.8 r c+3.4 +1.4 
'• y -3.8 -1.6 u -^-3.2 -1.1 " J - u -2.8 -0.9 

1B.2±|:S j£? Pb 11-6™ j*f pb 9.61U ±U Pb 



Table 3: First line: ti acceptance for the signal sample. Second and third line: Ntop and tt 
production cross section as a function of Mf 0p . The first error on each entry is the data 
statistical error; the second error on each entry is the sum in quadrature of all systematic 
errors listed in Table ||. 
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Figure 14: The combined Herwig (M top = 170 GeV/c 2 ) + VECBOS In(L*^°) distribution 
corresponding to the values of N top found (for M top = 170 GeV/c 2 ) in Table || (dashed) 
along with the data (solid). 
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6 Identification of b jets 

Each top quark event has two b jets. In contrast, direct W+ jet events contain b jets 
only at the level of a few percent |22[| . In this section we use two different methods to 
identify the b jets in the event (b tagging). In the first method, the Silicon Vertex Detector 
(SVX) is used to detect B hadrons by reconstruction of secondary vertices separated in the 
plane transverse to the beam from the primary interaction vertex as a result of the long 
B hadron lifetime. The algorithm used to reconstruct secondary vertices, "jet— vertexing" , 
and its performance are discussed extensively in Ref. O. For top quark events, the tagging 
efficiency (i.e the efficiency to tag at least one jet in an event as a 6-jet, including detector 
acceptance) is expected from Monte Carlo to be 24 ± 5 % in the signal sample and 19 ± 5 
% in the control sample. The efficiency is larger in the signal sample since more events 
fall within the fiducial acceptance of the SVX detector. From the number of top quark 
events derived from the ln(L*gJ°) shape analysis above, 1.5lj^g SVX tags are expected from 
top quark in the signal sample and 0.151q '^5 in the control sample. The expected number 
of SVX tags if no top quark were present in the sample is computed in the same way as 
in Ref. The dominant contribution is from Wbb, Wcc and mistags. The background 
estimate assumes that all the events are background and uses a tag probability for each jet, 
based on jet E^, r\ and track multiplicity, which is derived from a study of a large sample 
of inclusive jet data This contribution is found to be 0.47 db 0.06 events in the signal 
sample and 0.83 db 0.11 events in the control sample. Note that these estimates are derived 
directly from the data and do not rely on Monte Carlo predictions. Adding the other small 
background contributions to the tags (WW, WZ, Wc, Z —> tt and non-W events, see Q), 
the total expected number of tags assuming that the data contain no top quark events is 
0.58^o'o9 m the signal sample and 1.1 db 0.2 in the control sample. In the data, 4 events 
have a SVX tag in the signal sample (3 in common with the events selected in Ref. |TJ) 
and 1 event is tagged in the control sample. The probability that the tagging rate in the 
signal sample is consistent with the data being only background is about 0.4 %. On the 
other hand, the observed numbers of tags are consistent with the mixture of top quark and 
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background events expected from the kinematic analysis. Including mistags of top quark 
events and correcting the background estimates for top quark content in the sample, a total 
of 1.9 db 1.0 tags is expected in the signal sample and 1.2~^ ' 3 in the control sample. These 
numbers are summarized in Table 4, together with the probability that the observed tag 
rate is consistent with the respective expectation. 

As discussed in Ref. Q, b jets can also be identified by the presence of an electron 
or muon from semileptonic B decay. For this "Soft Lepton Tag" algorithm (SLT), the top 
quark tagging efficiencies are expected to be 19 db 3 % in the signal sample and 13 db 3 
% in the control sample. The expected number of SLT tags, assuming the data do not 
contain top quark, can be computed in a similar way as for the SVX jet-vertexing tagging 
algorithm. With all background contributions included, 1.2 d= 0.3 tags are expected in the 
signal sample and 1.4 db 0.3 tags in the control sample. In the data, there are 4 SLT tags in 



Sample 


Obs. 

tags 


Exp. tags 
backg. 


Exp. tags 
backg. + tt 


Prob. 
backg. 


Prob. 
backg. + tt 


Signal 


4 

SVX 


n r,o +0.12 
U - bS -0.09 


1.9 ± 1.0 


0.4% 


16.6% 


Control 


1 

SVX 


1.1 db 0.2 


1 2 +L3 

-0.3 


66% 


81.7% 


Signal 


4 
SLT 


1.2 ± 0.3 


2.2 ± 0.9 


4% 


20.7% 


Control 


1 

SLT 


1.4 db 0.3 


1 5 +L0 

L -° -0.3 


74% 


84.7% 



Table 4: Summary of b- tagging results in the signal and control samples. Also shown are the 
probabilities that the observed rate is consistent with background only, or with a mixture 
of top quark + background. 

the signal sample (3 of them in common with Ref Q ) and 1 (also in common with Ref. Q ) 
in the control sample. Two of the four events tagged by SLT in the signal sample are also 
tagged by the SVX jet-vertexing algorithm: one of these two is in common with Ref. (lj. 
In the signal sample there is again an excess of tags over the predicted background. 

The observed b-tags in the signal sample have a low probability of being entirely 
due to a fluctuation in tagging direct W+ jet background events. On the other hand, the 



35 



2 - 



C/3 

C 



> 

C 

CD 
> 

pq 



1.5 - 



1 - 



0.5 




Ln(L el ) 

Figure 15: Distribution in ln(L*gJ°) of the 4 events of the signal sample tagged by the SVX 
jet-vertexing algorithm. The expected tags in QCD VF+jet events (0.47 in total) are shown 
as a shaded histogram. 



hypothesis that the observed events are a mixture of background and top quark gives a 
good description of the observed tagging rates for both SVX and SLT tagging methods. 



7 Relative Likelihood of b— tagged and four jet Events 

The ln(L*gJ°) values of the signal sample events are listed in Table 5. One observes that 
5 out of 6 b-tagged events are at ln(L*^J°) > 0. The ln(L*^J°) distribution of the 4 SVX 
tags together with the dominant background from Wbb, Wcc and mistags, estimated from 
the inclusive jet parametrization, is shown in Figure 15. 

HERWIG predicts that about 80% of the top quark events (M top = 170 GeV/c 2 ) will 
exhibit a fourth jet in the CDF detector with transverse energy more than 15 GeV. Due to 
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SVX SLT 4th jet prim. lep. 
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3.1 
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43096-47223 
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42539-200087 
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• • * ii 


43351-266423 


1.1 
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45779-6523 
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* e 


42517-44047 
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• • fi 


44931-59686 


-4.3 


e 


47616-24577 


-5.0 
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Run-Event 


HKT) 


SVX SLT 4th jet prim, lep 


42913-59303 


2.2 


e 


45705-54765 


1.6 


• * e 


43276-101844 


0.2 


M 


45902-240098 


-2.2 


e 


46290-264893 


-2.9 


e 


45801-80320 


-3.4 


e 



Table 5: Summary of the signal sample events. The upper section lists events with at 
least one jet within the acceptance of the SVX tagging algorithm. The events in the lower 
section cannot be tagged by the SVX. The fifth column identifies those events which have 
four or more jets, and the last column identifies the lepton from the W decay as either an 
electron or muon. 

the small value of a s , the fraction of W+ > 4 jet events expected in a W+ > 3 jet sample 
is much less. Thus the requirement of a fourth jet should further enrich the sample in top 
quark events relative to QCD background. The signal sample contains 6 events with four 
or more jets with Ey(jet) > 15 GeV. These events are indicated with a * in Table 5. They 
are all at ln(L*^J°) > 0. Figure 16(a) shows their distribution in ln(L*gJ°) together with 
the prediction from the VECBOS W + 3 jets + HERPRT fragmentation routine. We recall 
that in this approach the fourth jet is produced by hard bremsstrahlung from initial and 
final state partons. Studies using W + 4 parton + SETPRT Vecbos simulated events yield 
similar predictions for background. The normalization chosen for VECBOS in Figure [Tc|(b) 
predicts 2.3 VECBOS events at ln(L*^°) < in Figure 16(a), while none is observed. This is 
compatible at the 10% C.L.. At ln(L*^°) > 0, 1.0 VECBOS events are predicted compared 
to the observation of 6 data events. The excess with respect to the QCD prediction at 
ln(L*gJ°) > already observed for W+ > 3 jet events in Figure 10(b), is therefore made 
relatively more pronounced by the requirement of a fourth jet, showing a positive correlation 
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between the ln(L*gj°) > and four-jet signature. The four-jet topology and b-tags are also 
strongly correlated. In three of the 6 four-jet events we find an SVX tag. The distribution 
in ln(L*gJ°) of these 3 events is shown in Figure 16(b). In the absence of top quark, 0.15 
background tags are predicted for the dominant direct W+ jet production, distributed as 
shown in the figure. 

A similar picture emerges for the SLT tag algorithm, since 3 SLT tags are identified 
in the W + 4 jet sample. In conclusion, 6 out of the 8 events at ln(L*^J°) > of Fig. 10(b) 
have at least one additional jet, and 5 of them are b-tagged. This is very unlikely to be 
due to background, and is more consistent with tt events. 

Using the methods described in Ref. (ij we have computed the top quark mass for 
the subset of events of the signal and control sample with exclusively four jets, by requiring 
Eta > 15 GeV and E^s < 10 GeV. Four events fulfill this requirement: three belong to the 
signal sample and one to the control sample. The three of the signal sample (all at m(L*i[ ) 
> 0) are in common with the W+ jet event sample of Ref. Q) and are among the 7 tagged 
events used in |l| for the derivation of the top quark mass. 

The masses of these 4 events are in the range 161 ± 11 GeV/c 2 to 172 ± 11 GeV/c 2 , 
lower on average but consistent with the result of Mt op = 174 ± 10 reported in Ref. Q. 
We have compared their mass distribution with the distributions expected for tt and direct 
W+4: jet events. The expected distributions for top quark and for direct W+ jet production 
are appreciably different. Within the very poor statistics, the distribution of the data events, 
shown in Figure 17, favours the tt hypothesis over QCD. 

8 Conclusions 

The kinematics of a sample of 49 W+ > 3 jet events was compared with the theoretical 
expectations for direct W^+jet production and tt quark pair production. It is determined 
whether a given W+ > 3 jet event fits better the expectations of direct W+ jet production 
as predicted by the VECBOS QCD Monte Carlo or top quark as predicted by the HERWIG 
Monte Carlo. The VECBOS predictions for W+ > 2 central jets and Z+ > 3 jet production 
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agree well with the observed data. A subsample of W+ > 3 jet events ("control sample") 
that should be enriched in direct W production events relative to top quark has been 
defined. VECBOS also gives a good description of the observed jet distributions for this 
sample. A separate subsample ("signal sample") is defined with the requirement that the 
three leading jets be central. It should be enriched in tt events relative to direct W+ jet 
events which form the main background. This signal sample contains 14 events. The jet Ey 
distributions for these events are unusually hard and not well described by the expectations 
from QCD and other backgrounds. By means of a suitable variable, ln(L*gJ°), events that 
are kinematically more top quark-like can be selected as those events with ln(L*]J ) > 0. We 
observe 8 such events, while we expect 1.7 from non-top quark processes. From a statistical 
analysis, which takes into account the systematic errors, we have derived a probability of 
0.8% for this excess to be due entirely to background fluctuations. The analysis was repeated 
for a number of different selection cuts defining the signal sample, and in the worst case a 
probability for such a fluctuation as large as 1.9% was found. A two component fit to the 
data that includes contributions from a 170 GeV/c 2 mass top quark and from QCD and 
other backgrounds gives a good description of the observed jet Ey distributions, and yields 
a tt production cross section of 11.6^g'7 j^io P^> consistent with the findings of Ref. 0. A 
similar two component fit to the background enriched control sample yields a cross section 
which is 1 sigma below this value, and statistically consistent with zero. 

With a secondary vertex 6-tag algorithm (SVX) we find evidence for bottom quark 
decay in four of the 14 events in the signal enriched sample. If the 14 events contained 
no contribution from top quark, only 0.58 events with such a secondary vertex 6-tag are 
expected. The probability for four events to be tagged due to a statistical fluctuation is 
0.4%. Similarly, this same event sample of 14 events contains 4 soft lepton tags (SLT) with 
an expected background of 1.2 events. The probability for four events to be tagged due to 
a statistical fluctuation is about 4% in this case. 

Additional information on the nature of the events at ln(L*gJ°) > was obtained from 
their large probability of containing a fourth jet. In the signal sample, out of a total of 8 
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events at ln(L*gT°) > 0, there are 6 four-jet events and 5 of them are b-tagged. Assuming 
that b-tags are indicative of ti pairs, one can argue that the events at ln(L*i[°) > show 
an increased top quark purity when the kinematic cuts are made more stringent (a 4th jet 
is required). We note that 5 out of 6 b-tagged events of the signal sample listed in Table 
5 are in common with the b-tagged sample of Ref. Q. This shows that, although the 
primary event sample selected in this analysis overlaps only in part with the W+jet sample 
of Ref. Q (25 events in common), the two analysis strategies have isolated the same physics 
process. The evidence for top quark reported in [|J was derived only on the basis of the 
observed excess of di-leptons and b-tags. The observation of a top quark-like component 
in the h^Ly™) distribution reported here provides additional evidence, independent of that 
provided by the counting experiments reported in Ref. Q , that our data contains a fraction 
of events more consistent with the decays of top quarks of mass around 170 GeV/c 2 than 
with the W+ jet background. 
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Figure 16: Distribution in ln(L*^[°) of W+ > 4 jet events, (a) data with VECBOS 
prediction, (b) events with a SVX secondary vertex and prediction (shaded) based on the 
observed jets in the events. 
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Figure 17: Distribution of the expected mass for VECBOS Monte Carlo events, analyzed 
as they were ti events (dashed histogram) (a) for the signal sample and (c) for the control 
sample cuts. Distribution of the preferred mass for tt (M top = 170 GeV/c 2 ) Monte Carlo 
events, analyzed as ti (b) for the signal sample and (d) for the control sample cuts. The 4 
events which allow the mass reconstruction are shown as solid histograms. 
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